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Since the discovery of half-filled cuprate to be a Mott insulator, the excitation spectra above 
the chemical potential for the unoccupied states has attracted many research attentions. There 
were many theoretical works using different numerical techniques to study this problem, but many 
have reached different conclusions. One of the reasons is the lack of very detailed high-resolution 
experimental results for the theories to be compared with. Recently, the scanning tunneling spec¬ 
troscopy (STS) [T] 12] on lightly doped Mott insulator with an antiferromagnetic (AFM) order found 
the presence of in-gap states with energy of order half an eV above the chemical potential. The 
measured spectral properties with doping are not quite consistent with earlier theoretical works. 
In this paper we perform a diagonalization method on top of the variational Monte Carlo (VMC) 
calculation to study the evolution of AFM Mott state with doped hole concentration in the Hub¬ 
bard model (HM). Our results found in-gap states that behave similarly with ones reported by STS. 
These in-gap states acquire a substantial amount of dynamical spectral weight transferred from the 
upper Hubbard band. The in-gap states move toward chemical potential with increasing spectral 
weight as doping increases. Our result also provides information about the energy scale of these 
in-gap states in relation with the coulomb coupling strength U. 


I. INTRODUCTION 

The spectral properties of the Mott insulators as a 
function of doping has been one of the key issues in 
studying the physics of high Tc cuprate superconduc¬ 
tors. There are many theoretical works on this topic in¬ 
clude exact diagonalization (ED)[3]-[5], quantum Monte 
Carlo (QMC)[6] method, dynamical cluster approxima¬ 
tion (DCA)[3, dynamical mean field theory (DMFT)[8]- 
[T2|, and also real-space Green’s function approach [13]. 
According to their results, generally, a clear spectral 
weight transferred from the upper Hubbard band (UHB) 
to lower Hubbard band (LHB) that situated at the chem¬ 
ical potential can be seen as doping increases. But the 
details are different. In Ref. laEHii], there are in-gap 
signals found. These signals become farther away from 
chemical potential together with the UHB as doping in¬ 
creases. There were also several results reported by using 
the cluster perturbation theory (CPT). For the electron- 
doped case [14], an in-gap state was found at the bottom 
of UHB, which is in agreement with the in-gap states 
seen by STS on Ca2Cu02Cl2[2] near an impurity. For 
the hole-doped case fTSHI^ . there were low-energy in-gap 
states with energy less than 0.2U. On the experimental 
side, the x-ray absorption spectra (XAS) [20l |2T] have 
observed the spectral weight transferred from UHB to 
LHB as doping increases but the broadness of the peaks 
makes it difficult to make detail comparison. 

Recently, STSP [2] reported for very underdoped 
cuprates with a long range AFM order probed the spec¬ 
tral function across the charge-transfer (CT) gap. It 
found some new results unexpected from earlier theo¬ 
retical works. In Ref. [T], finite density of states appear 
inside and throughout the CT gap and only a small en¬ 
ergy range at the chemical potential remains empty of 
spectral weight after holes are doped into the sample. 


The peak positions of these in-gap states at lower dop¬ 
ing can be above chemical potential up to 40% of CT 
gap. This is much larger than the in-gap states reported 
earlier [T4UT^ for HM calculations with U representing 
the CT gap. Besides the presence of in-gap states, when 
the system is doped with holes, there is also a system¬ 
atic evolution in the spectral weight distribution or local 
density of states (LDOS) measured at different positions. 
The in-gap states with larger spectral weight are situated 
more closer to the chemical potential, while at the same 
time the spectral weight of UHB moves to higher en¬ 
ergy. This relation is opposite to that found in earlier 
works HUHII], where in-gap peak moves to higher en¬ 
ergy as doping increases. Finally, the positions are anti¬ 
correlated between sites with higher spectral weight for 
UHB and in-gap states. This is consistent with an effec¬ 
tive doping picture. That is, at the position where the 
UHB has strong intensity, the effective doping is close to 
zero or no doping and the in-gap states are not seen. On 
the other hand, at the position where the in-gap states 
show up, doping is finite and the intensity of the UHB be¬ 
comes weaker. Note that similar results has been found 
previously in optical conductivity measurements [22]-[24], 
where there are also peaks around the scale of half an 
eV at low doping which moves to the lower energy as 
hole concentration increases. The discrepancy between 
these newly measured weight distribution and its doping 
dependence with earlier theoretical works has motivated 
us to examine the theoretical prediction again and more 
carefully. 

In this work we study the spectral evolution of Mott 
state with hole doping with a variational approach but 
with explicit presence of the AFM long range order as 
in the experiment [2]. We are particularly interested in 
the spectra of the unoccupied states above the chemi¬ 
cal potential. In the strong coupling regime {U > St), at 
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half-filling or in the parent compound, each site is already 
occupied by an electron with spin 1/2. Hence, when an 
electron is inserted into the half-filled state, it must cre¬ 
ate a doubly occupied site (doublon) and with final states 
in the UHB and there are no states inside the CT gap. 
However, after hole doping when an electron is inserted 
into the lattice, there are two possible final states. The 
original LHB splits into upper and lower spin density 
wave (SDW) bands by the presence of AFM order. In 
this case, finite in-gap spectral weight that corresponds 
to the upper SDW states shows up. These states, ac¬ 
cording to our calculation, behaves in a similar manner 
to the in-gap states recently found by STS[I1I1] with re¬ 
spect to the energy scale and evolution. These in-gap 
states also have components of states in the UHB with 
doublons despite the main contribution from upper SDW 
states. Thus these in-gap states are a mixture of LHB 
and UHB states and they absorb most of the spectral 
weight transferred from UHB. As hole doping increases, 
in-gap states move toward the chemical potential with 
increasing spectral weights and the energy separation be¬ 
tween UHB and LHB is effectively getting smaller. This 
provides a slightly different version from the ED result [3] 
without including the AFM order, which shows that the 
weights are transferred from UHB to LHB as the holes 
are doped but the band-edge separation has little depen¬ 
dence on the hole concentration. 

Below we first calculate the ground state of a one- 
band HM in the presence of the AFM order by means 
of the VMC method. Then several states with one elec¬ 
tron added to the ground state are proposed. These 
states contributing to the unoccupied states or the in¬ 
verse photoemission spectra (IPES) are orthogonalized to 
find the quasi-particle states. Then the spectral weights 
of these states are all calculated and compared with 
experiment [U [2] with respect to the energy evolution and 
the spectral weight redistribution. In addition, we also 
examine our results for different value of U to study the 
changes from weak to strong coupling. 


II. FORMALISM AND METHOD 

A well-known model which includes the low energy 
physics in Cu02 planes is the three-bands HM[25]. In 
this model, the parent compound without any extra 
doped holes has every Cu in configuration with a 
spin 1/2 hole. This is like a half-filled one-band HM with 
very large on-site Coulomb repulsion U and every site has 
a spin 1/2. When a hole is doped into the Cu02 plane, 
it resides at the Oxygen site [26]. Due to the strong su¬ 
per exchange interaction between the Cu spin and the 
doped hole on Oxygen, Zhang and Rice [27] found the in¬ 
teraction of two oxygen p orbitals and Cu d orbital lead¬ 
ing to three bands, the non-bonding states, anti-bonding 
triplet states, and the bonding singlet states known as 
the Zhang-Rice (ZR) singlet. There are large energy dif¬ 
ferences between the three states and only ZR singlet is 


assumed to be important for consideration. This ZR sin¬ 
glet in the three-bands model is similar to the vacant site 
when the hole is doped into the one-band HM. When the 
energy difference between ZR singlet and the Cu 3d^^ 
state, which is the effective CT gap, is not punitively 
large, the Cu hole can jump to its neighboring Oxygen 
to form ZR singlet while the original Cu turns into a 
3(i^^ configuration without spin. This is similar to the 
charge fluctuation process in one-band HM to turn the 
two nearest neighbor opposite spins into the short lived 
configuration of a doublon-hole pair. By making the cor¬ 
respondences of the doublon in one-band HM with the 
Cu 3d^^, the Hubbard gap with the CT gap and the va¬ 
cant site or hole with the ZR singlet, we could clarify 
the physics by studying the one-band HM instead of the 
more complicated three-bands mo del [28] [29]. 

The one-band HM we consider is 

H=-t y] {c\ aCj,a + h.c.) + y] (1) 

<i,j>,a i 

where t is the hopping integral of a single electron, < 
i^j > denotes the nearest-neighbour sites. U/t is usually 
taken to be 10 if not specified otherwise. 

The variational ground state we chose in the VMC 
method is the Jastrow type state with coexisting antifer¬ 
romagnetism and d-wave superconductivity [3Ql [3Tj . 

variational) = Pd—hPdl'^afm—ds ) ( 2 ) 

where Pd = is the Gutzwiller projection operator 
with d = di = representing the doublon 

number. The Gutzwiller factor g suppresses the dou¬ 
ble occupancy or doublon number when it is less than 
one. The Jastrow factor for doublon-hole binding [32] is 
Pd-h = njl - Qd-hdi rirp “ ^i+r)] where r connects 
the nearest neighbors and hi = (I — hip{l — hip is the 
number of holes on site i. This factor ensures the insulat¬ 
ing phase at half filling. This factor may come in different 
forms. Here we restricted the occurrence of free doublons 
that aren’t bound with holons with a variational param¬ 
eter Qd-h ^ 1- The wave function l^^afm-ds) with co¬ 
existing antiferromagnetizm and superconductivity has 
been proposed before [30]. 

\'^afm-ds)=P 

keMBZ 

(3) 

where restricts the state to have Ne electrons. The 
operators 

“L = «fc4.<T + 

4a = -^^fc4,a + «fe4+Q,a (4) 

correspond to the lower (a) and upper {/3) spin den¬ 
sity wave states (SDW) with coefficients = y(I — 
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Figure 1. The ground state staggered magnetization as a 
function of hole concentration for U=10t in a 12 by 12 lattice. 
The AFM order disappears around 0.18 doping in this model. 


—/==^), = 1 — My being a variational parameter 

proportional to staggered magnetization. Here we con¬ 
sider commensurate SDW and Q is chosen to be (7r,7r), 
and k is within the magnetic Brillonin zone (MBZ). 
The coherent coefficients Uk± and Vk± are defined by 




) and = 1 — , respec- 


.,2 = 1(1 _ 

^ y/A 2 +(£;fc±-/i )2 ^ 

tively. The plus/minus sign denotes upper/lower SDW 
states, Ek± = is the mean field SDW energy 

with Sk = —2t{coskx + cosky). The chemical potential 
/i is also taken to be a variational parameter. Finally, 
A/c = A(cos kx—cos ky) is the d-wave gap. For numerical 
convenience, our boundary condition is chosen to be pe¬ 
riodic in X direction and anti-periodic in y direction. The 
staggered magnetization obtained for the ground states 
is plotted as a function of hole concentration in Fig. 
The result that AFM order disappears around 18% agrees 
with Ref. [31] even though they have used a more sophis¬ 
ticated trial wave function [33]. By adding an electron to 
the ground state we can now calculate the IPES. We shall 
consider the simplest quasi-particle states and there are 
four kinds for each k point within the MBZ. 


|2fc,a) (6) 

i 

|3fc,a) = YE ~ Er n “ k+r+r)] 

i r r/—T 

(7) 

14,a) = “ k+T n “ k+r+r)] 

i r r/—r 

(8) 

where Iqn) denotes the ground state with particle num¬ 
ber N. States |1) and |2) both create an extra doublon 
in the ground state, so they belong to the UHB in the 
atomic limit. On the contrary, states |3) and |4) add 
an electron to a vacant site, and they belong to the 


LHB. Note that at low hole concentration, there are finite 
doublon-hole bound pairs generated by quantum ffuctu- 
ation. If we add an electron to the hole site bound with a 
doublon, it would create a free doublon. these states are 
also in the UHB, which we had confirmed by direct cal¬ 
culation of their energy. These states have large overlaps 
with |1) and |2) and they also contribute very little spec¬ 
tral weight which is proportional to the doublon number. 
So without loss of generality we shall exclude the process 
of creating free doublon from states |3) and |4) . 

To find the eigenstate within the chosen basis, for each 
k point in MBZ, we calculated the Hamiltonian ma¬ 
trix element by Monte Carlo algorithm {H{k^a))ij = 
{h,a\H\jk,a),{hj = 1^4). Since it is a non¬ 

orthonormal basis, we also need the metric tensor, 
{G{k,a))ij = {ik,a\jk,a)^ Next, we solve the 4 by 4 gen¬ 
eralized eigenvalue problem and obtain four eigenstates 
with an extra quasi-particle \'^)y_^^{k^a)) , (i = 1 ~ 4) 
and their energy cr). The energy to insert a 

quasi-particle is defined as 

ik,a) = £;]v+i {k,a)- EN+i,min (9) 

where the minimum eigenenergy EN-\-i,min is considered 
to be at the chemical potential. 

For each eigenstate at k, the spectral weight of insert¬ 
ing a particle contains two contributions Z+(/c,cr) and 
Zg(/c,cr), which are defined by 

z*+(fc,<T) = Kn+i(fc,4l4,j5iv)|' (10) 

Zi+{k,a) = ( 11 ) 


A similar procedure is also applied to study states with 
an electron removed from the ground state, this is for the 
photoemission spectra (PES). The states can be simply 
obtained by a transformation c. The basis states 

are 



(12) 


(13) 


i 


\^k,a) — (1 ^i+r+r)] 

i T r^—T 

(14) 

l^/c,cr) = ~ ~ (1 ~ ^i+r+r)] 

i r r^—T 

(15) 

After diagonalization we have four eigenstates with a par¬ 
ticle removed a)) , (i = 1 ~ 4) from the ground 

state and their corresponding energies are E]^_^{k, a) for 
each k. The energy to remove a particle becomes 

C~{k, a) = -Elf_^{k, a) + -E'jv-i.mm (16) 








4 



Figure 2. Spectral function for doping x = 0.0417,0.0833, 
and 0.125. A Lorentzian broadening with width F = 0.15t is 
applied to the delta functions for the eigenstates. The arrows 
indicate the value of as defined in the text. 


Similarly, the minimum eigenenergy is consid¬ 

ered to be at the chemical potential. The spectral weights 
to remove a particle are related to Z~{k, a) and {k,a) 
defined by 

Z^~{k,a) = \{^%_-^^{k,a)\c-k,cr\9N)f (17) 

Z^Q-{k,a) = \{¥j^_,{k,a)\c.k+QM9N)\^ ( 18 ) 

Finally we can combine the PES and IPES together, 

(k,cr)<uj-\-duj 

+ Z^+{k, a) + Z*- (fc, a) + zy (fc, a)) (19) 

The summation is over i, k and a. All results reported 
here are carried out on a 12 by 12 lattice. 


III. RESULTS 

The spectral function for U/t = 10 is plotted as a func¬ 
tion of energy for three hole concentrations in Eig. 

A sharp UHB peak around 9t is seen. More impor¬ 
tantly, a broad band of in-gap states appears in the range 
0.17/ ^ 0.357/ which we believed to be mostly unoccu¬ 
pied upper SDW states. These states are about the same 
energy range as the newly found signals in the STS[T1|2]. 
As hole concentration increases, the UHB weight shifts 
toward higher energy while the in-gap states moves to¬ 
ward the chemical potential, i.e. the energies of in-gap 
states decrease. The spectral weights of these in-gap 
states also increase with doping. There is clearly a spec¬ 
tral weight transferred from UHB to low energy states 
(LES) that are between chemical potential and UHB. To 
verify the relationship between in-gap states and the up¬ 
per SDW states, we calculate the inner product between 
these states. The upper SDW states can be constructed 
in the same way as wave functions |3) and |4) except 
now we restrict the electron to be inserted into the up¬ 
per SDW band defined in Eq. 0- 



Figure 3. of the upper SDW states as a function of mag¬ 
netization. The energy together with magnetization are col¬ 
lected from four different doping x = 0.0417,0.0833,0.125, 
and 0.1667, here magnetization is inversely proportional to 
hole concentration. 

li’u-SDwikja)) = 

Hi n <1 - (20) 

r T^—T 

Where {i\Pk)(j = —crhke^^^^ + At each 

of the k points, we found large overlap (> 0.8) between 
\'^u-SDw{k^<j)) and the in-gap states. Thus these in¬ 
gap states are essentially the upper SDW states although 
there are contributions from states \ lk,a) and |2/c^cr) which 
are in the UHB. Besides these states near the chemical 
potential, there are also contributions from lower SDW 
states that are now vacant due to hole doping. The en¬ 
ergy scale of these states is roughly determined by the 
coupling between the states \3k^a) and |4/c^cr)- Consider¬ 
ing the transition between upper and lower SDW states, 
our result gives a possible explanation of the half an eV 
peak in optical conductivity measurements [22H24] that 
shows a decreasing absorption energy and an increasing 
weight with more doping. This will be left for future 
works. 

In the presence of AFM long range order, staggered 
magnetization opens a gap between the upper and lower 
SDW states. In our case, the effective staggered magne¬ 
tization is proportional to the variational magnetic field 
My. To illustrate this relation, we define by the lowest 
eigenenergy of the quasi-particle states that has an inner 
product with I'lpu-SDwik^ <j)) larger than 0.8[34], which 
provides a good indicator of the lower edge of the in-gap 
states. Thus, positive correlation is expected between 
the upper SDW band edge and AFM strength (M) as 
shown in Fig. Since (M) is inversely correlated with 
doping, this gives a natural explanation of the reduction 
of energies of in-gap states as doping increases. Next 
we shall examine the spectral weight transferred from 
UHB to LES. We calculated the total weight Wuhb and 
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Figure 4. Total spectral weight for UHB and LES. In under¬ 
doped regime, Wuhb and Wles evolve linearly. By fitting, 
we found the slope to be —1.66 for Wuhb and 2.69 for Wles- 

Wles- 

Wuhb = ^ Z^+ {k,cj) + Z^+ {k, a) 

^*+(A:,cr)>4t 

Wles= Z^+{k,a) + Z^+{k,a) ( 21 ) 

^*+(/c,cr)<4t 

Wuhb and Wles different doping are shown in Fig.|^ 
This agrees quantitatively with the previous results from 
ED [3], ED + cluster DMFT [35], and CPT[T8] despite the 
fact that AFM order was not considered in these works. 
This shows that the appearance of AFM order doesn’t 
affect total spectral weight transferred as the sum rule 
should be satisfied. 

If we consider each site at the atomic limit, UHB is 
completely unoccupied at half-filling, hence it has a spec¬ 
tral weight equal to 1. This weight Wuhb reduces to 
1 — X when x holes are doped into the system. Upon 
doping, an electron could be added to the empty or hole 
site in two choices from the spin degrees of freedom, 
hence the spectral weight for Wles is 2x. However, it 
is known that beyond this atomic limit there should be 
a dynamical correction [36| that comes from the coupling 
between these states which enhances the weight transfer 
and would give Wuhb = I — x — a and Wles = 2x -\- a. 
According to our calculation, the renormalization a at 
U/t = 10 is around 0.67x. 

Since the AFM order depends on U, it is important to 
examine the evolution of spectral functions with different 
values of U. In Fig.|^a), the spectral functions for insert¬ 
ing a quasi-particle in the ground state are plotted for a 
range of U values. For U/t > 6, the separation between 
UHB and LES is clear and as expected, UHB energy 
scales with U. The trend suggests that as U becomes 
smaller, weight of in-gap state at x = 0.0556 becomes 
larger and their component of |1) and |2) increases. For 
weak or intermediate U/t the weight of the in-gap state is 
comparable or even larger than that of UHB. To further 
examine the U dependence of these spectra, we plot the 
upper SDW band edge as a function of U in Fig. [^b). 
Before U reaches Uc ^ St to enter the Mott region, the 
^0 is proportional to U/t. This is expected from a mean- 
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U/t 


2.0 
^ 1-5 
^ 1.0 
0.5 

^ 0%0 0.01 0.02 0.03 0.04 0.05 

<M>/(U/t) 

Figure 5. (a)IPES for different U at x = 0.0556. {h)^o/t at 
different U. {c)^o/t as a function of {M)/{U/t) for U > 8t. 
This shows that the virtual exchange interaction J accounts 
for the results. 


field treatment for a weak or intermediate coupling U in 
the one-band HM, as the gap due to AFM order is pro¬ 
portional to {M)U. Once the system enters the Mott 
region with U greater than Uc, it is the superexchange 
interaction J that determines the AFM order. Hence 
^0 ^ {M)J ^ {MW{U/t). ^0 plotted as a function of 
{M)/{U/t) in Fig. [^c) for U > Uc- This shows that one 
could expect a maximum energy for in-gap states when 
U/t is near the critical value for a Mott transition. 

Now lets examine the spectral weight for LES. As 
shown in Fig. [^a), Wles is extrapolated to zero at 
half-filling for U much larger than Uc- It has a sudden 
increase below U/t = 6. This suggests that these in¬ 
gap weight won’t disappear at half-filling in intermediate 
coupling case. Moreover, as U decreases to the weak cou¬ 
pling regime, as shown in the U/t = 4 case in Fig. ia), 
the major peak appears at the location of upper SDW 
states. Those upper SDW states begin to merge with 
the UHB and the excitation gap becomes an SDW gap 
and there is no clear separation between UHB and up¬ 
per SDW states. This is in accordance with the mean 
field theory. Therefore, the disappearance of SDW states 
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(a) 



X 



Figure 6. {d^)WLES for different U at x — 0.0417, 0.0833, and 
0.125. (b)Slope oiWhES for different U taken from data at 
X — 0.0417 and 0.0833. 


at half-filling as U increases is a characteristic of Mott 
transition. Furthermore we can also examine the slope 
of the WleS‘ An important feature of Mott insulator is 
the spectral weight oi 2x a for the LES. In the infi¬ 
nite U limit, according to Ref. [36], a is proportional to 
1/U. In this regime, as U decreases, there is an increase 
in the slope of Wles- fhe free particle limit, on the 
other hand, there is no UHB and all of the weight 1 x 
is distributed close to chemical potential. The transition 
between strong and weak coupling can be clearly seen in 
Fig.gb)., starting from 7//t = 16, as U decreases, there 
is an increase in slope oiWLES- As soon as U crosses the 
Mott transition (I/c), the slope has a sudden drop. At 
7//t = 6, the slope is much less than 2, suggesting that it 
has crossed Uc and Mott physics no longer applies. 


IV. CONCLUSION 

In summary, by using a variational approach to study 
the HM in the AFM phase, we construct several quasi¬ 
particle states to study the evolution of PES and, in par¬ 


ticular, IPES, with hole doping. The substantial amount 
of spectral weight inside the Mott gap is due to the mostly 
unoccupied upper SDW bands. These states interact 
strongly with the UHB so that there is a large spectral 
weigh transferred from UHB to these in-gap states. Al¬ 
though we use one-band HM instead of specific oxygen 
orbitals, we are able to capture the detail evolution of 
the spectrum with doping observed in the recent STS on 
cuprates. Our results also agree with previous numeri¬ 
cal works with respect to the dynamical spectral-weight 
transfer with a renormalization of about a ^ 0.65x at 
U/t = 10. This agreement is a bit surprise as previous 
works have not included AEM order [3l [35]. Maybe this 
is due to the fact that for IPES, there is a much smaller 
incoherent spectral weight [37] and quasi-particle states 
we considered have almost all the spectral weights. 

We like to point out that in the lightly doped cuprate 
samples [T] in experiments, there is no density of states at 
the chemical potential. Since the sample is quite inho¬ 
mogeneous, we believe that there is a strong localization 
effect that depletes the density of states at chemical po¬ 
tential. Since we have not considered the localization or 
disorder in our calculation, our result has a small peak 
near chemical potential when an electron is added to the 
lower-SDW states. Eortunately, the in-gap states that we 
focused on has energy much larger than this localization 
gap, so that we can account for them without consider¬ 
ing the disorder. An interesting possibility of localization 
may be due to the checkerboard pattern observed in[T]. 
In the strong coupling limit of the HM or the t — J model, 
it is recently shown [38] that at very low doping, there are 
states with checkerboard patterns involving SDW. Hence 
our account of SDW as the source of in-gap states may 
be a reasonable approximation. Eurther work on this is 
underway. 

Our result also shows that the spectral distribution 
has a non-monotonic behavior when U is increased above 
Uc to enter the Mott physics. At half-filling, there is a 
small SDW gap at small U, and also a finite spectral 
weight at low energy (Wees) above chemical potential. 
But the gap changes to the much larger Hubbard gap 
as U becomes larger than Uc and there is absolutely no 
spectral weight within the gap. After doping the energy 
of in-gap states changes from mean field type (oc U) to 
t-J (oc 1/U) type and left a peak around Uc- At small 
doping the low energy spectral weight above chemical 
potential is proportional to doping for U > Uc- These 
properties might be useful to observe the Mott transition 
by applying pressure. 
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